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MM Lala TT HH HH HH HH $$ II NN NNNN HH HH 
MM MM TT HH HH HH HH $$ II NN NN HH HH cece 
Lal Lal TT HH HH HH HH SS II NN NN HH HH eves 
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LL I] SS 
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; H floating Hyperbolic Sine routine 


HISTORY ; Detailed Current Edit Wi story 
D CLARATIONS ; Declarative Part of Module 
MTHSHSINH = Standard H floating SINH 
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Q: Sf ayes Macro V04-00 Page 
2 CMTHRTL.SRCIMTHHSINH.MAR; 1 


1:4 
1:2 

~TITLE MTHSHSINH : cating Hyperbolic Sine routine 

IDENT /1-006/ : File: MTHHSINH.MAR EDIT: RNH1006 


SREARAAEAEAAAAAAAEAAAAAAAAAAAAAAAAAAARAAAAAAAAAAAAAAAAAAAAAAAAAAEAAAAEAEAAEEAEEE 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL to eM _GORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESER 


THIS SOFTWARE IS A ght UNDER A LICENSE AND MAY BE USED AND gt 
ONLY IN accom NCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
anol ya ton OF THE ABOVE pith ol NOTICE. THIS SOF TW ARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


*® ® 
*® ® 
*® ® 
*® * 
*® * 
*® a 
*® ® 
*® ® 
*® ® 
® ae 
s® TRANSFERRED. * 
** ® 
® ® 
** * 
*® ca 
*® ® 
** * 
:® ® 
: a 
*® x 
** ® 
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einusion ; H floating Hyperbolic Sine routine ib a 14 -1 1384 
H 
( 


3 


——_ —©& 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eeepokat ite NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


LAR RRSAAAAAAALALALLALALALALALALALALE ALE AAA RARE SASS RRA R AAA RRL AAA RRA A ARAL | 


s FACILITY: MATH LIBRARY 
; "ABSTRACT: 
: MTHSHSINH is a function which returns the H floating hyperbolic sine 


; of its H floating point argument. The call is standard 
; ; call-by-reference. 


SOSOSoOOoOoCoOoOOSooOSoOoOoOoOoOSoOSoSoOoOSoO 


SOOCOSOOOOOSOOSOOSOOSOSSOSOOS OOOO OOOOOoOoOSoSoO 


; VERSION: 1 
; HISTORY: 

cies A. Wheeler, 11-Sep-1979: Version 1 
: ; MODIFIED BY: 


DONAUVE WN $$ O OD NAUE WIN OOD NAUE WN @§ OOONAUE WN OWOONOUS Un 


PERLE LES FWWWWIWIWIIINIUINONononononononononn 3 2 3 oe 


SOoOSSSCSCSOOOSOSOSOSOSOSOOOSOSSOSOSSSOSSOSSOSOSOSOSOSoOSOSoOSoOOoOSOSOSoO 


SOOCCOCOOOOOCSOOSOSOOSOSOSOOSOOOOOOOCOOoOOSoOoOoOO 


SoOooooooooooooo 


Loat ing Wyperbet fc Sine reunites” = See 138% 1:30:07 AX/VMS Macro V04 
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; Detailed Current Edit History EP-1984 MTHRTL. SRE RCIM MTHHSINH.MAR; 1 


-SBTTL HISTORY ; Detailed Current Edit History 


; ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
; Edit History for Version 1 of MTHSHSINH 


Adapted res MTHSGSINH vere tee 1-001. JAW 11-Sep-1979. 
Use MTHSHEXP_R6,. SBL 3-O0ct-1979 

Make H_0. ; a iterel. Set “B-0¢t-1979 

Don't store reserved operand before signal. SBL 7-Feb-1980 
changed lower Limit for Chebyshev approximation from 2**-56 to 


ftintines ageane call to EXP for input values between 57*Ln2 
Changed att t'pel floating point divisions by 2 to interger 
subrtracts of 1 from the exponent field. 

Extended maximum range to 163584*Ln 

cnaneged « Lo ye for computing EXP(!xi=ln2) to reduce error. 


Changed Shsinh to @hsinh(AP) in error logic RNH 22-Sept-81 


SO te 


COOOCCCOCGCOOOOOOOOOCOOOOoO 
NOUS WN OOONOUE Wt" 

ee ey 

@eees 

Cooooo 

ooooo 

USWwn— 

eeeuene 


ow 


COOQOQooooooooooooooooo 


SOOoOooooooooooooooooooo 
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1-006 beCLART Hone {*per bolle Sing routine aut Gosepoione Sissies PATANTe eeesarenciqu.mans: 2% ¢% 


§ -SBTTL DECLARATIONS 3; Declarative Part of Module 


5 : INCLUDE FILES: 
> EXTERNAL SYMBOLS: 

“EXTRN MTHSSSIGNA 

“EXTRN. MTHSK FLOOVEMAT 
SEXTRN MTHSHEXP_R6 


EQUATED SYMBOLS: 


+ MACROS: none 


+ PSECT DECLARATIONS: 


COQOCOOQCoQooQoooooooooooo°o°oo 


SOOOOCoCoOOCOOOCOOCOoOoOOoOOooOoooo 


FP PLMWAMWMVININ | “GDOOOOOCCOSOSOOSOSOOSOSOSOSOSOOSOSOOSOOOOOOOOOoOO 


i 
i 
7 
8 
0 
¢ 
4 
5 
9 j 
0 8 
0 89 
0 90 
0 91 
i 
00000000 «9 .PSECT _MTHSCODE PIC, SHR,LONG, EXE ,NOWRT 
+ a 3 program section for math routines 
0000 97 OWN STORAGE: none 
00 9 98 ; 
00 99 ; 
0000 100 ; CONSTANTS: 
0000 101; 
0000 10¢ 
0000 103 H_57_LO0G_2: 
6DB93AB1 3€134006 99 104 [QUAD *X6DB93AB13C134006 
C70997E6 FF9D90B9 00 105 [QUAD *XC70997E6FF9D90B9 : 57*ln2 
99 106 H_16383_L0G_2: 
ESEQAGSE 62DE400E 0010 107 QUAD “XE 3EQAGSEG2DE4O0E 
68251812 €66964DB 0018 108 SQUAD *X6B251812C66964DB : 16383*Ln2 
wsoearer ezeccgne Bh 199 "-165-L055 sapearerezee 
07666780 Liars) ° iW Saeed Qua Re aeea eters t : 16384*Ln2 
A39E2FEF $2€44000 0 118 rae “A396 26 EF 6644000 na 
40006730 93C7F 0 114 cies sQUAD *%4000673093C7F 35 : (high 99 bits of in2)+2"#-99 
FO CODOBF 9D o 116 7 7” “.QUAD “XFOCB9SOBCODOBF9D 
ors 3559 D674CD98 $f 9 QUAD oNSP SeO029D674CD98 ; Un2 = H_LOG_2_HI 
09 119 HSINHTAB: 
c377E97C 425456 AD 3 9 -LONG X42 3436 AD OXC377E97C s 6.506912239471978471647518358407046E -26 
536 FAES 7elest ee 6 ; -LONG “Hol ¢3r 88 6. R. HB cee AES : 3.867997522833381902976185781667616E -23 
BCAAt 303 Lent ? 4 -LONG a Bet :  1.957294395551957859061134358006348E -20 
al 48s grascobs $ -LONG uct cB, teas :  8.220635243500749521323757431341128E -18 
0€3077 aa 0 128 LONG X950C3 bb. *x0C3D77 :  2.811457254347797192926287843860811E -15 


4 
(3 
3 
-8 
“6 
=4 
-3 
-1 
0 
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7.647163731819805091328739180656633E -1 
1.605904383682161463834581885720363E -10 
2.505210838544171877496282899225692E 
2.755731922398589065255745053025947E 
1.984126984 126984126984 126972636750E 

8 .333333333333333333333333333385565E 

1 .666666666666666666666666666666657E 

0 .000000000000000000000000000000000E 
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Standard H floating SINH 
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Brperses ts Sine seuttos yp oe =13 4 01:40:07 VAX/VMS M 


EP 

EP-1984 11:25:42 (CMTHRTL.S 
-SBTTL MTHSHSINH = Standard H floating SINH 
+ 
FUNCTIONAL DESCRIPTION: 
HSINH = H floating point function 
HSINH{X) is computed as: 


If 1X! < 2-57, HSINH(X) = X, 

IF 2ee-57 =< {X: <0.25, HSINH(X) = Ghebyehev $ 
If 0.25 =< Ix! < are ts rh ge eF = (HEXP(X) = 

If 5?#ln2 =< IX! < 1638%#Ln2, HSINH(X) = si n(x 
If 16383*Ln2 =< 1x! < 16384*(n2, then HSINH(X) 

If 16384*ln2 =< iX!, then ouertine. 


CALLING SEQUENCE: 
hsinh.wh.v = MTHSHSiNH(x.rh.r) 
-or- 
CALL MTHSHSINH(hsinh.wh.r, x.rh.r) 


04-00 
RCIM THHSINH.MAR; 1 


Because an H-floating result cannot be expressed in 64 bits, it is 


returned as the first argument, with the input 
to the second argument, in accordance with the 
Standard. 

INPUT PARAMETERS: 


LONG = ; define Longw 
x = 2 * LONG ; Contents o 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 

hsinh = 1 * LONG 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 


; SIDE EFFECTS: 


; Signal: MTH$_FLOOVEMAT if leyesetnd . Xi ig # reserved yc tess i RO/R3 
. (copied to the signal mechanism vector CHFSL 2 Margt 4 abt 


a ee ee ee me ae me a a ee ee De a me ee ee a ed a ed ad dd od dd 
ROPOPONoNoONOnY 2 2 
SOOOOOOO OO OOOO OOO 0000 60 09 09 09 09 09 09 09 SI NIN SIS NSN IPA AA AAA AAOA MMIII 
AALAND $9 OD NAUES WIN 9 OO NAME WIN  O OONAU EWN OC OD NAU EWN 0 ODNOAOUE WN 


; Associated message is: ‘FLOATING OVERFLOW IN “RATA LIB 


arameter displaced 


ord multiplier 
x is the argument 


wag Mk is reserved 


; operand -0.0 unless a user supplied (or any) error odes “changes CHFSL_MCH_RO/R1 


; NOTE: This procedure disables floating point underflow, enables integer 
; overflow. 


04-00 Page 6 


MTHSHSINH ro_V 
1-006 MTHRTL.SRCJMTHHSINH.MAR; 1 (4) 
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ing Hyperbolic Sine routine 16-SEP-1984 01:40:07 VAX/VMS Macro 
: Standard H floating SINH aes 91 358:03 l 

1 

1 

1 


ENTRY MTHSHSINH, “M<IV, R2, R3, R4, RS, Rb, R7, RB, RY, R10> 

: standard call~by-reference entry 

3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 
3; set handler address to jacket 
; handler 


Nm—OOOon 


6D O0000000'GF 9 


3; case of an error in routine 

MOVAB G*MTHSSJACKET_HND, (FP) ; Set handler address to jacket handler 
: If an error, convert signal to user PC 

; and resigna 

> RO/R3 = ix 


6D GQ0000000'GF 9€E 


50 _08 BC 70FD MOVH ax(AP), RO : X! = @value(AP) 
57 50 7DFD MOVO RO, R7 
50 8000 8F AA BICW2 #*x8000, RO : RO/R3 = 1X! 
SFFF BF 30 B1 CMPW 0, #*xSrFF + compare {X! with 0.25 
23 «18 BGEQ  GEQ_T0_0.25 : branch if !X! >= 0.25 
: IX! < 0.25 
3FC8 BF 50 —s«B1 CMPW RO, #°X3FC8 zs compare {Xi with 2**-57 
06 18 BGEQ  GEQ_TO_2m57 : branch if {X!} >= 2*#=57 
IK! < 208-57 
04 BC 57 7DFD MOVO R7, @hsinh(AP) ; Result = argument 
04 RET 3; return 


> 208-57 =< 1X! < 0.25 


NNAAAAAAAA AO MUUIUUUIUUIUUIUIUIUIE BE EE EERE EEE SWWIWIoonororonororononononnns «= Zo 


OOM 600060 09090909 S99 SF SF APNINPININPININNIN SO OOOOE PMN OUD OO OOOOONMNINNINMOCOCO Ie 


PUPP BS BS BS BS BB BE BE ENA NNW NINN PONNONPONONPUNININ 2 sO 
CONOAOME WR OOO NOU EWN OOO NOAUE WN 9 OONAUES WN" OOONOULSW 


Porononononononononononononononononononononononononononononononononononornononoronononononony 


1 : 
1 
1 GEQ_TO_2m57: 

50 50 64FD 01 MULH2 0,R0 ; Get ARG**2 for POLYG 

FEF3 CF OC 50 75FD 01 POLYH RO, #HSINHLEN-1, HSINHTAB 

1 : RO/R3 = SUM(Ci*X**i), with 
1 ; last coefficient zero 

50 57 64FD 1 MULH R7, RO : MULH2 | ARG, and then 

048C 50 57 61FD 01 ADDHS R7, RO, @hsinh(AP) : add in ARG with overhang and return 
04 ! RET ; with result in first argument 
1 3 
; 3 0.25 =< ix! 
1 a 
FEA2 CF 50 71FD 0 6£6_10.0,635 RO, H_16383 G2 compare {Xi with 16383*ln2 
e m tAr 
4514 a BGTR OTR. TRAN, 16385 COG_2 3 branch if ixi > 16883eine 

1 3 
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MTHSHSINH ; H floating Hyperbolic Sine vention 16-SEP-1984 01:40:07 VAX/VMS Macro V04-00 Page 7 
1-006 MTHSHSINH = Stendard H floating SINH ae tt 9 91:39:08 YRTHRTL. SREIMTHHSINH.MAR: 1 ° (4) 
17 59 ; 0.25 =< IX! =< 16383*Ln2 
1 get 
FEBA CF 30 71FD Q17 66 CMPH RO, H_57_L0G_2 ; Compare iXi to 57*ln2. If greater 
14 f o BGTR ONLY_ONE_ TERM ; only one call to HEXP is necessary. 
17 65 ; 0.25 =< IK! < 57#Ln2 
17 08 3 
17 6 
50 57 72FD 017 68 MNEGH R7, RO ; RO/R3 = - 
n° sal F 16 017 69 JSB MTHSHEXP_R6 ; RO/RS = HEXP(-Xx) 
f 7DFD i 70 MOVO RO, -(SPY 3 p sh HEXP(-X) on stack 
5 7DFD 1 71 MOVO P é O/R3 =X. 
00000000 "EF 16 018A Le JSB MTHSHEXP_R6 : RO/R3 = HEXP(X) 
50 E 62FD 0190 7 SUBH P)+, RO 3 RO/R3 = HEXP(X) = HEXP(-x) 
50 01 A2 0194 27% SUBW2  #*x0001, RO + RO/R3 = (HEXP(X)<HEXP(=X))/2 
04 BC 50 7DFD 0197 75 MOVO RO, @hsinh(AP) : Return SINH(x) in fisrt argument 
04 19C 6 RET j 
B10 77; 
19D 78 ; S7*ln2 =< IX! < 16383*Ln2 
019D 79 ; 
019D 80 
019D 281 ONLY_ONE_TERM: 
00000000 EF 16 019D Hf JSB MTHSHEXP_R6 3 RO/R3 = HEXP(iX!) 
57 53FD A 8 TSTG 7 ; Check sign of X 
04 14 O1A6 84 GTR POSITIVE ; If negative change sign of 
50 50 S2FD O1A8 85 MNEGG 0, RO : HEXPCTX!) 
O1AC 86 POSITIVE: 
50 01 A2 OQ1AC 87 SUBW2 #*x0001, RO : RO/R3 = sign(X) *HEXPCiX!)/2 
04 BC 50 7DFD O1AF 288 OVO RO, ahsinh(AP) : Return SINH(x) in fisrt argument 
04 0184 89 RET ; 
0185 90 ; 
8182 91 ; 16383*in2 =< ix! 
1B5 8 ; 
1B5 9 
1B5 294 GTR_THAN_16383_L0G_2: 
FE65 CF 50 71FD 0185 95 CMPH RO, H_16384_L0G_2 ; Compare iXi_to 16384*Ln2 
23 18 0188 96 BGEQ ERROR s if 16384*ln2 =<ixXi, overflow occurs 
018D 97 
018D 98 
01BD 99 ; 
01BD 00 ; 16383*ln2 < iki < 16384*Ln2 
4 01 ; 
50 FE6E CF 62FD 188 88 SUBH2 H_LOG 2HI, RO : RO/R3=:Xi-H_LOG 2 HI | 
00000000" EF 16 OIC 04 JSB MTHSHEXP_R6 3 RO/R3 = HEXPC({XT-A_LOG_2_HI1) 
7 73FD O1C9 05 TSTH R7 3; test the sign of X 
2 18 1C€C 8 BGEQ 10$ ; branch if X >= 0 
50, 50 726 Oc MNEGH RO, RO RO/RE=s ign X) *HEXPC:X; -H LOG_2_HI1) 
4 50 FE69 CF 65FD 01D be 10$: MULH H_LOG_ 2_L0, RO, R4 > R4/R3=HEXP( {Xi -H_LOG_2_HT) *H-LOG_2_L0 
04 8C 50 54 63FD 01D 0 SUBH RZ, RO, ahsinh(AP) : Return HEXP(iX!-Tn2) in first 
04 ODF 10 RET 3 argument 
1e0 11 
10 $1¢ | 
1E 1 
1E 14; 
1E 15 ; 16384*ln2 =< iXi, error 
| 
| 
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16 ; 
1 
18 ERROR: MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; cengieten value 

ASHQ #15, #7, R RO/R3 = result = poserves operand -0.0 
4 
8 


MTHSHS INH : 
1-006 artisn 


50 1 OF 


rc so error handler can mod 
00000000" EF ef FB 
8? BC ror? 


cLRO P y the result 
CALLS 4 MTHSSSIGNAL signal error and use real user's PC 
RO, a@hsinh(AP) Restore result 


; return 


| 


ee a a a a ad ——- 
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MTHSHS INH 3; H floating Hyperbolic Sine ro ste 16-SEP-1984 01:40:0 
Symbol table “ eure green 1obe ST ids 


ERROR 1€0 R 1 
GEQ -10 35 g3 168 R 1 
GEQ"TO- 152 R 1 
GTR-THAN_16383_L0G_2 1B5 R 1 
HSIAH = 4 
HSINHLEN s 
SINHTAB R 1 
H_16383_LO0G_ 10 R 1 
eae ee 
Nt oe-$-to 9000040 RO 
( = $o000006 
MTHSSJACKET_HND eeeeeeee =X 01 
MTHS$SIGNAL eeneneee =X 8600 
MTHSHEXP_R6 eereeeee§=6X 8600 
MTHSHSINA 00000120 RG 01 
MTHSK_FLOOVEMAT eeeeeeee = X 00 
ONLY_GNE_TERM 00000190 R 4 
ve QOOOOIAC R 1 
= 0000000 
$ ewer ener omen eece a 
! ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
- ABS. 00000000 9.3 0.) NOPIC USR CON ABS LCL 
_MTHSCODE OO0001F7 ¢ 50 1.) PIC USR CON REL LCL SHR_- EXE 
aenaneeeecasecceessonsen} 
H ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 30 00:00:00.0 0:00:01 .38 
sounens processing 114 BBR! 0. Bs hoes 
s 1 86 3:8 :01. 0:00:03.54 
Syabol table sort 0 0:00:00.01 0:00:0 3 
Pass 2 71 00:00: 9.78 8:8 702.85 
Symbol table output 4 00:00:00. :00:00. 
Psect synopsis output 3 8 :00:00.01 0:00:00. 
Cross-reference output :00:00. 8 0:00:00. 
Assembler run totals 309 00:00:02.6 0:00:12.4 


The working set Limit was 900 pages 

5105 bytes (10 pages) of virtual memory were used . buffer the intermediate code. 

There were 10 pages of symbol table space ot lecatee & hold 21 non-local and 7 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to “define 1 macro. 


7 VAX/VMS Macro v04-00 Page 
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—— NOEXE NORD NOWRT NOVEC BYTE 
RD NOWRT NOVEC LONG 


ee en ‘" 
| 
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Sry esis toe 
VAX-11 Macro Run Statistics 


fowmetan ne nernraen enema ean an cece $ 


! Macro Library statistics ! 


tw moamannrene rene se see eoee ee 


Macro Library name Macros defined 
_$255SDUA28: (SYSLIBJSTARLET.MLB;2 ereamaacat 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHSINH/OBJ=OBJ$:MTHHSINH MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 
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